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ABSTRACT 

We develop a numerical formalism for calculating the distribution with energy of the (internal) pairs formed in a 
relativistic source from unscattered MeV-TeV photons. For GRB afterglows, this formalism is more suitable if the 
relativistic reverse-shock that energizes the ejecta is the source of the GeV photons. The number of pairs formed 
is set by the source GeV output (calculated from the Fermi-LAT fluence), the unknown source Lorentz factor, and 
the unmeasured peak energy of the LAT spectral component. We show synchrotron and inverse-Compton light- 
curves expected from pairs formed in the shocked medium and identify some criteria for testing a pair origin of 
GRB optical counterparts. Pairs formed in bright LAT afterglows with a Lorentz factor in the few hundreds may 
produce bright optical counterparts (R < 10) lasting for up to one hundred seconds. The number of internal pairs 
formed from unscattered seed photons decreases very strongly with the source Lorentz factor, thus bright GRB 
optical counterparts cannot arise from internal pairs if the afterglow Lorentz factor is above several hundreds. 

Subject headings: methods: numerical - radiation mechanisms: non-thermal - relativistic processes - shock 
waves - (stars:) gamma-ray burst: general - (stars:) gamma-ray burst: individual (GRB 
130427A) 


1. INTRODUCTION 

The LAT instrument onboard the Fermi satellite has detected 
a high-energy emission at 100 MeV-100 GeV, extending well 
after the GBM prompt phase, for dozens of GRB afterglows 
(first Fermi catalog - Ackermann et al 2013). The properties 
of the LAT afterglow emission are: fluence above 100 MeV 
$ = 10 _5±1 erg cm -2 , light-curve peak at 10-20 s after trig¬ 
ger, post-peak flux decay vF v oc f _1 ' 3±0 3 monitored up to 1 
ks (sometimes longer), photon spectrum C v oc u~ 2J ±0 - 2 . 

The isotropic energetic output of the brightest LAT after¬ 
glows, Ery = 10 53±1 erg, is 10-100 percent of the GRB output 
at < 1 MeV, the corresponding number of afterglow photons 
above e 7 ~ 100 MeV being N 1 = E 1 /e 1 ~ 10 57±1 . The 
fraction f± of these photons that form pairs depends strongly 
on the Lorentz factor T of the medium that produced the LAT 
afterglow emission, because T determines the lab-frame colli- 
mation of photons and the threshold energy for pair-formation, 
and on the source radius R, which sets the optical-thickness to 
photon-photon absorption. Taking into account that R ~ V 2 ct, 
with t the observer time, it follows that f± has a very strong 
dependence on T. 

From the escape of the higher-energy LAT photons (~ 10 
GeV), Abdo et al (2009) have set lower limits T 0 > 200 — 1000 
for the Lorentz factor of several LAT sources during the prompt 
emission phase (burst). Consistent with that, Panaitescu et al 
(2014) found that, for F > 200, photon-photon attenuation 
does not yield a spectral signature but, for V < 200, the at¬ 
tenuation of photons above 1 GeV should be detectable. Ad¬ 
ditionally, as shown in this article, for T < 500, the number 
of pairs N± = f±N 1 is higher than the number of electrons 
energized by the forward shock and, for T < 250, is larger 
than the number of ejecta electrons energized by the reverse 
shock. Therefore, the formation of pairs from LAT photons is 


of importance at least for those GRBs/afterglows whose LAT 
spectrum displays a spectral softening at GeV energies. 

Assuming a single power-law for the LAT spectral compo¬ 
nent, Panaitescu & Vestrand (2014) have calculated analyti¬ 
cally the emission from the pairs formed only from photons to 
which the GeV front is optically-thick, leading to the conclu¬ 
sion that pairs can account for the brightest optical counterparts 
(flashes) observed during the prompt phase. A larger number of 
pairs, but of lower energy, are formed by the photons for which 
the LAT emission is optically-thin (to pair-formation). In this 
work, we calculate numerically the distribution with energy of 
all pairs formed from high-energy photons (assuming a broken 
power-law spectrum), integrate it over the deceleration of the 
blast-wave that produces the LAT afterglow, and track numeri¬ 
cally the pair radiative cooling, to obtain accurate light-curves 
for the optical flashes produced by pairs formed in GRBs and 
afterglows. 

The effects arising from scattering of the high-energy pho¬ 
tons on (cold) electrons existing in the source or on the already- 
formed pairs are ignored. Such scattering increases the source- 
frame photon escape path, which increases the probability that 
any photon forms a pair and the total number of pairs formed. 
As we shall see, for sources that are optically-thin to photon 
scattering, scattering on the already-formed pairs occurs with a 
smaller probability than pair-formation thus, to a good approx¬ 
imation, the effect of photon scattering on the pair-formation 
rate can be ignored. The same is true for scattering on the ejecta 
electrons, if the Lorentz factor of the GeV source is in the few 
hundreds. However, if that Lorentz factor is or exceeds sev¬ 
eral hundreds, scattering on ejecta electrons should occur with 
a higher probability than pair-formation. In this case, the num¬ 
ber of pairs formed from unscattered seed photons, calculated 
below, underestimates the true number of pairs. 

We also ignore the formation of (external) pairs ahead of 
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the afterglow blast-wave, which loads with leptons the ambi¬ 
ent medium and accelerates it (Beloborodov 2002), changing 
the afterglow dynamics (Kumar & Panaitescu 2004). Conse¬ 
quently, the formalism presented here for the emission from 
pairs is more pertaining to a GeV source that is located well be¬ 
hind the forward-shock, so that most pairs form in the shocked 
fluid and not ahead of the blast-wave. That condition points to 
a relativistic reverse-shock as the origin of the LAT afterglow 
emission (as could be the case for GRB 130427A - Panaitescu 
et al 2013). 


2. PAIR-FORMATION IN A RELATIVISTIC SOURCE 


2.1. High-energy spectral component 


The number of pairs formed at any observer-frame time t 
over a dynamical timescale is derived from the observable 0.1- 
10 GeV fluence <1», the spectrum of the high-energy afterglow 
emission (with the 0.1-10 GeV spectral slope being the only ob¬ 
servational constraint), the source redshift z, and the unknown 
Lorentz factor T of the high-energy source. At the redshift of 
the afterglow source, the afterglow photon spectrum is assumed 
to be a broken power-law: 


dNj _ dN 1 ( (e/tb) a e < £b 
de de eb \ (e/eb)~ 13 £b < e 


(1) 


where ft, is the spectral-break energy at redshift z (i.e. the 
peak of the luminosity vL v spectrum), a and u being the 
low and high-energy spectral slopes. For synchrotron and 
inverse-Compton emissions, a has four possible values: 2/3 for 
optically-thin synchrotron (sy) or inverse-Compton (ic) emis¬ 
sion from un-cooled electrons, 3/2 for optically-thin sy/ic from 
cooled electrons (i.e. with a radiative cooling timescale shorter 
than the age of the source), -1 for self-absorbed synchrotron 
emission (but is unlikely that the source magnetic field is suf¬ 
ficiently high for self-absorption to be important at MeV), 
and 0 for the inverse-Compton scattering of self-absorbed syn¬ 
chrotron emission. Then, the spectral slope around 2 measured 
by LAT above 100 MeV indicates that £b < 100 MeV (in the 
observer frame), jd ~ 2, and that the low-energy spectrum is 
not observed, being dimmer at 10 keV-10 MeV than the GRB 
spectrum. 

The normalization factor of equation (1) is simply set by the 
measured fluence $ 
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5.10 52 (z + l) 2 cm is the luminosity distance and 



2.2. Peak/break energy of the LAT component 

A lower limit on the observer-frame break-energy £b can 
be set by requiring that the 0.1-10 keV afterglow emission 


measured during the X-ray light-curve plateau (at 0.3-10 ks) 
by Swift/XRT, of about T xr t — 10 -10 ergcm _ 2 s _1 (O’Brien 
et al 2006), is not dimmer than the extrapolation T\ a t of the 
GeV afterglow spectrum, for an afterglow of GeV fluence 
$ = 10 ~ 5 ergcm ~ 2 at f = 10 s, decreasing as $ oc t ~° 3 , 
and with a high-energy slope /3 = 1.1: 

= 3.10- 8 (jL) 11 ergon-V 

(4) 

For a = 2/3, the high-energy spectral component does not 
overshine the X-ray plateau flux if ej,(l ks) 50 keV. For 
a = 3/2, the high-energy spectral component is dimmer than 
the X-ray plateau if £b(lks) > 10 MeV, but could be the 
X-ray plateau if £b — 1(//1 fcs) _1 ' ±0 -' MeV. This evolu¬ 
tion is consistent with the f -3 / 2 expected for the peak energy 
of the forward-shock synchrotron spectrum, but implies that 
£b(10 s) > 1 GeV during the burst, which is inconsistent with 
LAT observations, that do not show a high-energy component 
peaking in the LAT window. 

However, X-ray plateau measurements are often lacking dur¬ 
ing the early GeV afterglows monitored by LAT, thus the above 
low limits on £5 cannot be derived for individual afterglows. 

The sub-MeV burst light-curve may also set a constraint 
on the unknown et, in the following way. For St, ~ 10 
MeV and a = 2/3, the LAT spectral component yields 
a 100 keV flux T lat = (<S>/t grb )(0.1 MeV/e b ) 2 ~ a = 

2 . 10 ~ 9 $_ 5 t~ 3 b y£b 7 ^ 3 er § cm _ 2 s _1 during a t gr b = 10 s burst 
(using the notation X n = X(cgs)/ 10” and measuring photon 
energies in eV). This emission is sufficiently below the typ¬ 
ical flux of a bright burst, T gr b = 10 ~ 5 ergcm~ 2 s _1 , that 
it does not overshine a fast-decaying F gr b oc f - ( 2 -=-4) tail 
(O’Brien et al 2006). In contrast, an energy-break £b that 
falls below 100 keV will produce a burst emission Ti a t = 
($/tgrb){OAMeV/£ b ) 2 ~ 0 = 10 - 6 $_ 5 f“ 1 6 il ergcm- 2 s - 1 

(independent of for a high-energy LAT spectral slope fi ~ 2) 
that rivals that of the prompt emission. 

Thus, a bright LAT afterglow following a slowly-fading GRB 
may have a break-energy £b < 100 keV, but one following a 
burst with a steep decay must satisfy £& 2S> 1 MeV during the 
burst tail. 


2.3. Optical thickness to photon-photon absorption 


For an isotropic distribution of photons in the frame of the 
shocked fluid, which moves at Lorentz factor I in the lab-frame 
(at redshift z ), the optical thickness to a photon of energy e! a is 


T 77 ( £ o) 


0~ e 

47ri? 2 



dN^f^e’^') (5) 


where primed quantities are in the shock-frame, cr e is the Thom¬ 
son cross-section for electron scattering, 

R(t) ~ (z + l)~ 1 ctT 2 (6) 


is the source (shock) radius at observer-time t (corresponding 
to the arrival-time of photons emitted by the visible edge of the 
source; equality in equation above holds for an undecelerated 
source; a factor 4/3 applies to the right-hand side for a blast- 
wave decelerated by a wind-like medium), 9' is the angle of 




Fig. 1.— Optical-thickness to pair-formation as function of observer-frame photon energy for a source of high-energy photons with the indicated parameters 
(fluence in cgs units) and for three values of the break energy eb . 


incidence between the test-photon of energy e' a and a target- 
photon e', d,P/dO' = (1/2) sin0' is the probability of two pho¬ 
tons interacting at an angle 9' 


4(40') = 


2 (m e c 2 ) 2 


(1 — cos 9')e' a 

is the threshold-energy for pair-formation, and 

/ 77 (4'> 9') = ^ (2 + 2x ~ 2 - x~ A ) In (a; + \J x 2 - l) 
— (l + x ~ 2 ) \/l — x ~ 2 
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is the cross-section for photon-photon absorption, with 


which leads to an integral over 9' that can be calculated eas¬ 
ily. Dropping the integral over 9' and assuming that most pairs 
form at threshold are further simplifications that lead to the cor¬ 
rect dependence of r 77 on source parameters: 
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where 


e t h(e) = 4r 
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is the lab-frame threshold-energy for a lab-frame incidence an¬ 
gle 9 = r _1 -C 1 (the source motion at T collimates photons 
within an angle T" 1 around the direction of motion). Using 
equations (1), (2), and (2), one arrives at 


'1 e'e' 
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2 m 2 c 4 
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The integral in equation (5) is calculated numerically; Figure 
1 shows r 77 (e) for a photon of observer-frame energy 


£ = 


Fe' 


z + 1 z + 1 


( 10 ) 


corresponding to the typical relativistic boost (F) of a photon of 
shock-frame energy e'. To extract the dependence of r 77 (e) on 
the source parameters I’, $, z and observer time t, an approxi¬ 
mation to equation (5) is needed. Extending the approximation 
/ 77 ( x) = (3/4) ln(2cc)/a; 2 , accurate for x 1, to all x > 1, 
allows the second integral in equation (5) to be calculated ana¬ 
lytically, but the resulting integral over the incidence angle 9' is 
not so nice. The dependence of r 77 on source parameters can 
be obtained by setting / 77 (a;) = const and by approximating 
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(14) 

having switched to observer-frame photon energies (equation 
10). Using equation (7), the conditions above become 


£ b < £th (£) 
£th{£) < £b 


£ < £ 
£ < £ 


withe = [4T 2 (?n e c 2 ) 2 ]/[(z + 1)4] = 4.6 Z 2 T 2 3 £ b l GeV, 
where Z = (z + l)/3. 

Equation (14) shows the obvious fact that r 77 is proportional 
to the photon column density (N^/R 2 oc 4>/f 2 ), and that it has 
a strong dependence on the source Lorentz factor. Given that 
a > 0 and /3 > 2, equation (14) also shows that r 77 (e) in¬ 
creases with photon energy for e < e and decreases with it for 
i < e, with the maximal optical thickness reached at e. 

The coefficients missing in equation (14) depend on the pho¬ 
ton spectrum slopes a and (3. Figure 1 shows optical thickness 
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for the representative values a = 2/3 and /3 — 2, obtained nu¬ 
merically by integrating equation (5). In the asymptotic power- 
law regimes, the numerical approximation is 


'77 


00 = < 


°- 57 d £ 9 

1 2.3^1 


£■ < £ 


:.3 16 

The two branches above intersect at 


O c Z 10 / 3 $_5 _ —1/3 = _ 

r 10/3 ,2 4/3 £ 11 £<£ 

2.3 ^ 1^6 


(15) 


Thus, perfect power-law LAT spectra set only a lower limit 
on the source Lorentz factor: T > 300 t t 1 / 3 . Such a weak 
dependence on the afterglow fluence $ and epoch t of obser¬ 
vations suggests that the measurement of curvature in the LAT 
spectrum would yield a fairly accurate determination of T. Ob¬ 
viously, the non-detection of the high-energy afterglow emis¬ 
sion is not necessarily proof of high absorption and does not 
tell us anything about T. 

2.4. Total number of pairs 


e = 10 z~ 2 ^ GeV 

£ b,6 


(16) 


The total number of pairs is an integral over the photon spec¬ 
trum of the absorbed fraction g[r 77 (e) 


where the optical thickness is maximal: 


r 77 (e) = 5.7 


£ 4 $-5 
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Equation (15) can now be written as 


T ni £ ) ~ T ~n( £ ) 
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From here, it follows that 
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For /? = 2, the optical-thickness to photon-photon absorption is 
independent of e b for photons of energy lower than e (i.e. pho¬ 
tons with threshold energy above e b ), hence the lower limit e_ 
above which the photon-front is optically-thick is also indepen¬ 
dent of e b , as illustrated in Figure 1. 


N ± = 



de 
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de 
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with the photon spectrum of equation (1). To calculate the frac¬ 
tion of absorbed photons corresponding to the optical-thickness 
to pair-formation r 77 (e) (equation 15), consider a medium of 
geometrical thickness A and linear absorption coefficient a, 
and in which the production and absorption of photons is ho¬ 
mogeneous (same at any location).Then, the fraction of photons 
that are absorbed is 


9=[ ^(l-e"“< A ->) = 

JCI A T 


r/2 r«l 

1 t » 1 


(24) 

obtained by integrating the photon absorption from the medium 
inner edge (x = 0) to its outer boundary (x = A), and with 


r = aA. 


In the case of pair production in a decelerating source, the 
photons radial distribution is not uniform. In this case, the frac¬ 
tion of absorbed photons is 


9 = 



(25) 


where dn-y/dy is radial distribution of photons normalized by 

/g 1 dy(dn-y/dy) = 1 and 


r (y 


1) = f a{z)dz = t(0 — > 1) f ^j^dz (26) 
J v J y dz 


Equation f 19) can be reinterpreted as 

r 77 (e)<l if r > r± = 310 1/2 ^ 5 /4 (22) 

For T < T±, the photon-front is optically-thick to photons in 
the (e_, e+) range, which widens with decreasing T (see equa¬ 
tion 21), covering the entire LAT window if e_ < 100 MeV, 
which is equivalent to T ^ Ti at = 125 Z^ 1 _/ 3 t 1 1 ^ 3 . For 
such low Lorentz factors, the LAT emission is heavily absorbed 
and the afterglow high-energy emission undetectable. At the 
other extreme, if T ^ 2.2 Vi„t, then e_ // 10 GeV and the 
LAT emission is weakly absorbed. For T; ai < T < 2.2 T; at , 
the LAT emission is moderately absorbed, photon-photon ab¬ 
sorption rendering a spectrum that curves downward at higher 
energies, for a power-law intrinsic spectrum (e.g. figure 2 of 
Panaitescu et al 2014). 


is the absorption optical thickness from coordinate y = x/ A to 
the outer edge at y = 1, r(0 —> 1) = r being the entire optical 
thickness of the medium. Substituting in equation (25), we get 



Therefore, as long as the absorption coefficient a is propor¬ 
tional to the density of the to-be-absorbed photons, the fraction 
of absorbed photons from a medium depends only on the opti¬ 
cal thickness r = J a(z)dz of that medium and does not ’’care” 
about the exact spatial variation of a. 
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The integral in equation (23) is calculated numerically; for an 
analytical estimate, we use the approximation given in equation 
(24). 

Case 1. For £b > ib/2 (equation 19), the maximal optical 
thickness (equation 17) satisfies r 77 (e) < 2, hence r 77 (e) < 2 
for any photon. With two branches for the photon spectrum 
(equation 1) and two for the optical thickness (equation 14), the 
integral in equation (23) splits in three integrals 


N ± = 


dN-y 




for the more likely case £& < e. Using equation (17), this con¬ 
dition requires that £b < 100 Z _ 1 r 2.3 MeV, which is satisfied 
by LAT spectra and which, together with the working condition 
£b > ib/ 2 , requires that l > 1 where 


r = 


100 z- 


$ 


1/5 

-5 


2/5 


t 


(29) 


The scaling of the integrals in equation (28) with e is a . 
In e, and e~^~ a \ respectively; taking into account that (3 > a, 
this implies that most pairs are formed from (the second inte¬ 
gral, corresponding to) photons with energy above the spectral 
break £{, > £i,/2 ~ few MeV and below the energy for maximal 
optical thickness £ ^ 10 GeV, interacting with photons above 
threshold energies of about 1 GeV and 1 MeV, respectively. 

For a = 2/3 and /3 = 2, equation (28) yields 


N± 2e = ib 10 55 f 0.67 + 0.22 In - 5 ^-^ Z 8 -|^- (30) 

V -*£ 6,7 / 2.3 

Using equation ( 6 ), the optical-thickness to photon scattering 
of the pairs is 

$ 2 — 5 
4pl0 
l 1 2.3 

(31) 

Ignoring the logarithmic term, this implies that the pairs are 
optically thin (to photon scattering in the Thomson regime, be¬ 
cause most pairs are cold) for T > F T with 

F t = 138 Z^tf 2 ' 5 (32) 
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(33) 




for the more likely case £b < £- (for F > F, this is implied by 
the working condition £ b < £b/2). 

The integrals in equation (33) show that most pairs form from 
(the second integral, corresponding to) photons above the spec¬ 
tral break £& < ib/2 ~ 3 MeV and below £_ ~ 1 GeV, for 
which the photon front is optically thin, interacting with pho¬ 
tons above threshold energies > 1 GeV and few MeV, respec¬ 
tively. 

For a = 2/3 and (3 = 2, equation (33) leads to 


iV ± “*=<“ 10* (l.l + 0.11 m J|4L) (34) 

After calculating the pair optical thickness to photon scattering 
as done above for £& > ib/ 2 , it can be shown that the minimal 
Lorentz factor for optical-thinness is close to that in equation 
(32). 

Equations (30) and (34) show that the number of pairs 
formed has a weak dependence on the (unknown) break £/, of 
the photon spectrum, varies like + 2 (as expected for a two- 
photon interaction), and has a strong dependence on the source 
Lorentz factor, resulting in part from the dependence of the 
threshold energy for pair formation on F and in part from the 
decrease of the photon density with source radius (which is pro¬ 
portional to r 2 ). 


2.5. Scattering of afterglow photons on internal leptons 

It is worth comparing the number of pairs with that of elec¬ 
trons existing in the two possible source of GeV afterglow pho¬ 
tons, the forward and reverse shocks. 

For a Wolf-Rayet GRB progenitor with a mass-loss rate 
dM/dt = 10 -5 Mq /yr, blowing a wind of terminal velocity 
v w = 10 3 knr/s, the wind baryon density is 


_ dM/dt 
( ^ 4t Tm p v w R 2 


3.0 x 10 35 
R 2 


cm 


(35) 


The ratio of the number of formed leptons TV; = 2TV± (eqs 30 
and 34) to the electrons energized by the forward-shock is 


Case 2. For £b < £/,/2, t 77 (e) is relative to 2 as in equation 

( 20 ) but with £_ larger by a factor 2 1 /( /3-1 ) than j n equation 

(21) and £ + smaller by a factor 2 1 ^ 1_a ^ than in (21). Having 
two branches for the photon spectrum and three for the optical 
thickness, the integral of equation (23) splits in four; 


TV± 


ciTV 7 min 



dN-y 

de 
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Sb 



TV; 

Nfs 


N± 

irnR 3 


~ 2600 


t?r|. 3 


(36) 


taking into account that the wind-like medium is made of 
elements heavier than hydrogen (with one electron for two 
baryons). Thus, for V < Tf s = 540 the pairs 

are more numerous than the forward-shock electrons. 

If we assume that i) the reverse and forward shock baryons 
contain about the same (kinetic plus thermal) energy and ii) 
the ejecta are normal matter (with one electron per baryon), 
then the number of ejecta electrons is at most T times larger 
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than that of the forward-shock’s (as in the case of a semi- 
relativistic reverse-shock): N rs < r Nf a . Then, the number 
of pairs exceeds that of the ejecta electrons for T < 1 ’, ., = 
270 Z$% 9 f" 1/3 . 

The above suggest that emission from pairs is of importance 
for GeV afterglow sources with a Lorentz factor in the few 
hundreds, but pairs may radiate at a different energy than the 
reverse or forward-shock electrons, where the pairs could dom¬ 
inate the afterglow emission even if they are fewer. 

It is also worth investigating if scattering of pair-forming 
photons on existing (reverse and forward-shock) electrons or on 
the already-formed pairs could change significantly the number 
of pairs formed from unscattered photons. Most photon scat¬ 
tering occurs on leptons that are cold. That is certainly the case 
for the pairs, most of which are born cold (see the distribu¬ 
tion of formed pairs with energy in Figure 2, left panel), and is 
likely true for the (reverse-shock) ejecta electrons and the ambi¬ 
ent medium electrons (swept-up by the forward-shock) because 
they should be cooling fast radiatively if their synchrotron emis¬ 
sion were to account for the observed GeV afterglow. 

For a scattering optical thickness r sc , the effective photon- 
photon attenuation thickness is r 77 = -\/r 77 (T 77 + r sc ), there¬ 
fore scattering on cold leptons is negligible when r 77 (e) > 
t sc . As shown in §2.4, most pairs form from photons with 
energy e < i, for which r 77 (e) is that in the first branch 
of equation (15). The optical thickness to Thomson scat¬ 
tering on reverse-shock electrons (which are more numerous 
than in the forward-shock) is t sc ^ a e (Y Nf s ) / (Att R 2 ) ^ 
5.10 _ 3 iJ/(r 2 3 ^ 1 ). Thus, r 77 (e) > t sc is satisfied for e 
£rs = 10(r 2 . 3 /^) 5 (^/^- 5 )MeV. 

The optical thickness t± to Thomson scattering on (already 
formed) pairs is that given in equation (31), thus r 77 (£) > T± 

is satisfied for£ > £^ )7 = 50(Z/r 2 .3) 4 ($_5/ff) MeV. Scat¬ 
tering on cold electrons of observer-frame photons with en¬ 
ergy above Ek n = Tm e c 2 /(z + 1) = 34 Z - 1 r 2 .3 MeV oc¬ 
curs in the Klein-Nishina regime. Given that c^r <; £f~ n , it is 
worth considering scattering on pairs in the KN regime, when 
the scattering cross-section akn(s) — (3/ 8 )cr e ln( 2 x)/x with 
x = e/£kn- In this case, the condition t 77 (£) > r±(e) is satis¬ 
fied for £ > e£$ = 25(Z/r 2 . 3 )3/ 2 ($i / 5 2 Ai) MeV. 

To the above identification of the photon energies for which 
scattering increases significantly the photon escape path and at¬ 
tenuation, we add that, according to equations (28) and (33), 
most pairs are formed from photons with energy in the range 
(£;,,£) (for£{, > £;,/2) or (£;,,£_) (for£{, < if,/ 2). This means 
that scattering is important for photons in the lower part of those 
energy intervals and not so important in the upper part. Equa¬ 
tions (28) and (33) also show that each decade of photon energy 
provides an equal contribution to the number of pairs. There¬ 
fore, scattering has a negligible effect on the number of pairs if 
the logarithmic length of the upper interval is larger than that 
of the lower interval. Using the expressions for i, £_, e rs , and 

£pair, it can be shown that scattering on pairs should not in¬ 
crease much the number of pairs if T > 120 (for £5 > if,) 
and r > 140 (for £/, < if,), while scattering on reverse-shock 
electrons is negligible if T < 370 (for £ 1 , > if,) and T < 780 
(for Sf, < if,), having left out sub-unity powers of the param¬ 
eters Z, $, t. Adding to these that the reverse-shock electrons 
are more numerous than the pairs if F > T rs = 270, it fol¬ 


lows that scattering of the pair-forming photons below the LAT 
range increases the total number of pairs only through scatter¬ 
ing by reverse-shock electrons and only if T is at least several 
hundreds. 


2.6. Pair distribution with energy 

The shock-frame energy of a pair e' v depends on the ener¬ 
gies of the incident photons, test-photon e' 0 and target-photon 
e', the incidence angle 9 ', and the center-of-momentum (CoM) 
frame angle (j>" at which the electron and positron emerge, 
measured relative to the direction of motion of the photons’ 
CoM. In the shocked-fluid frame, the CoM moves at velocity 
P'cm = (e'o + e') /(e'o + e'), the corresponding Lorentz factor 
being 

(’ _j_ p 

T' cm = ; 0 = (37) 

\/2e' 0 e'(l — cos 6’) 

In the CoM frame, the incident photons have the same energy 


e"= - e ' o e'(l-cos (0 


Ik(CO') 


collide head-on, and form an electron and a positron of equal 
energy e", moving in opposite directions, at angles (j>" and 
7 r — (j)" relative to f3' cm . Then, the shock-frame electron and 
positron energies are 

e' ± = r' cm (e" ±p"cP' cm cos <t>") (39) 

where p"c = e" 2 — m 2 c 4 is the electron/positron momentum 

in the CoM frame. 

To obtain the distribution of formed leptons with their shock- 
frame energy, dNf/de'±, we start from equation (23) in shock- 
frame photon energy e' a 

dNf __ o dN , _ dN-t f r 77 (e^) r 77 (e^) < 1 

dZ ~ 2 ^Z 9M£ ° )! - i 2 r 77 (e')»l 


where the factor 2 accounts for two leptons being created from 
one photon. The term r 77 ( f ' 0 ) of the first branch offers a way 
to calculate dNf/de ±, by expanding it (as in equation 5), lead¬ 
ing to the differential pair-number in the 4-dimensional space 

K,e',e',r] 


d 4 N[ 

de'dQ'de'dQ" 


cr e dN 7 dP dN 1 dP t fJ j a ,^ 
AttR 2 de’ dCV de' dW 111 ^ 1 } 


where dP/dil" = sin0"/2 (because pairs emerge isotropi¬ 
cally in the CoM frame), with / 77 the photon-photon absorp¬ 
tion cross-section of equation ( 8 ). For an isotropic distribution 
of the incident photons ( dP/dVL' = sin0'/2) 

_ d' Ni _ cr e dNj , dNj , , , 

de' 0 dd'de'd(cos (j)") I&ttR 2 de' Q de' ' 77 0 

(42) 

From equation (39), de'± = ±T' cm f3' crn p"c(d cos ej>"), thus 

d 4 Nf _ d 4 Nf dcoscj)" 

de’ 0 d9’de'de' ± de' o d0'de'd( cos (j)") de'± 


I 67 tR 2 de' 


dN 7 gin q, 44 M e 'o e '' 6 ') 


de' r cm P' cm p"c 
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Then, the distribution of leptons with energy is the integral of 
the differential pair-number above over the energies of the inci¬ 
dent photons and over the incidence angle 


dNi 

de'± 


167ri ? 2 


de 


, dN. 


/•IT 


de' 


7 f dO 1 smO' 
o 



de 


,dN-y 

sfW 


f~n(e' 0 e' ,0') 


TOgC 4 ) 


(44) 


where the limits e' min (e±) and e' max (e'±) on the integral over 
the spectrum of target photons are determined from 


r' 


the shock-frame photon energy corresponding to the observer- 
frame £. A 7 m e c 2 < e! b branch does not exist because Co < 4 
photons form pairs in interaction with photons of energy e' > e' b 
and the corresponding pair energy e± ~ e '/2 is in one of the 
branches above. The above distribution was derived for e' b < 
(which requires T > f) but is also correct for e' b > e' (requiring 
that T > T) with e b and e! swapped. 

For £b < if,, we have t 77 (£) > 1 for £_<£<£+ and 

t 77 (£) < 1 otherwise, thus 



7 1 d b < jm e c 2 < e!_ 

7 _/3 e’_ < 7 m e c 2 < e' + 

7 -(fl+i-a) 7 m e c 2 


(48) 


together with e' > e'thiCn Unfortunately, the term 
F'cmftcmP" leads to a fourth degree equation in e' that cannot 
be solved analytically to obtain the integral limits e' min (e'±) and 

^max( e ±)‘ 

Those limits can be calculated numerically and used to inte¬ 
grate equation (44), with the following two corrections. First, a 
multiplicative factor g[r 77 (ep)]/ T 11 {e' 0 ) should be applied to 
the e' Q integrand, to account for the correct absorption frac¬ 
tion g{£' 0 ). That ensures that no more than dN 1 {e' 0 ) pho¬ 
tons form pairs. Second, a multiplicative factor min{l,iV 7 [> 
e tti( e o)]/[-W 7 (> e l 0 )g(T 11 {e' 0 ))\} ensures that the number of ab¬ 
sorbed test photons iV 7 (> e' 0 )g[T^(e' 0 )] does not exceed the 
number of target photons iV 7 [> e' th (e' 0 )] above the threshold 
for pair-formation. 

The slopes of the pair distribution with energy dNi / de' ± that 
results from integrating equation (42) can be inferred if the 
crude approximation e± ~ e' 0 /2 is made. This approxima¬ 
tion is suggested by that pairs emerge most likely at an CoM- 
frame angle <f>" = 7r/2, hence e’± = T' cm e" = (e' a + e')/2 
(from equation 39), and by that most pairs are formed from a 
test-photon of energy e' a larger than the e' of the target-photon. 
The latter is suggested by that most pairs are formed in interac¬ 
tions with target-photons close to (but above) the threshold for 
pair-formation (the integrand in equation 5 shows that r 77 oc 
e , (d/V 7 /<ie , )/ 77 (e , ) oc (e') 1_ “ ln( 2 a;)/a ; 2 with x 2 oc e' a e\ thus 
r 77 oc (e') _Q with a > 0 ) and by that e' Q e' th {e' 0 ) (optical 
thickness t 77 (£) is maximal at photon energy e £5 = £th(s) 
- equation 17). 

The approximation e' ± ~ e' 0 /2 implies that dNi/de'± oc 
dN,y /de '±. Then, equation (40) leads to 


dNi dNry , dN.y . f r i t 

~de ^ 9 ^ T77 ( e± )^ ^ mm! 1 , T 77 (e ± )| (46) 


with an approximation for the absorption factor g that has the 
correct dependence on r 77 . From here, the slopes of dNi/de' ± 
can be easily calculated using the photon spectrum (equation 1 ) 
and the optical thickness (equation 14). 

For ib < £b, we have t 77 (£) < 1 for any photon, thus 


dNi ( 7 1 e' h < 7 m e c 2 < e! 

(I7 * [ 7-(/ 3 +t-a) g/ < ^ meC 2 


(47) 


with e' the shock-frame photon energy for which r 77 is max¬ 
imal (given in equation 17), 7 = e' ± /(m e c 2 ) the pair’s ran¬ 
dom Lorentz factor in the shock-frame and e' = e(z + 1)/F 


for the more likely case e' b < e'_, i.e. for F > f (equation 29). 
If 0.63 f < T < f, then e'_ < e b < e' + and the first branch 
above is 7 “° with e b and e'_ swapped. For T < 0.63 f, we 
have < e 'b and, in addition to the preceding case, the second 
branch above is 7 -1 , with e' b and eV swapped. 

Equations (47) and (48) indicate that the pair distribution 
with energy has up to three power-law branches, with four 
possible values (— 1, — a, — /3, — /3 — 1 + a) for the slope of 
each branch. Figure 2 shows that the pair distribution ob¬ 
tained numerically by integrating equation (42) displays only 
the highest-energy branch for £{, < £/, because the range of en¬ 
ergies over which the lowest-energy branch occurs is too nar¬ 
row, from 7 = 1 to 


7(e / ) 


0 + !)g = 15 T2.3 

2rm e c 2 Z£b ,7 


(49) 


For £b <C £;,, all three branches are found in the numerical pair- 
distribution. The lowest-energy branch is short, extending from 
7 ~ few to 


7(0 = 


(z + l)e_ 


= 26- 


r 5 t 2 

1 2.3 6 1 


2rm e c 2 Z 5 $_ 5 ' 

while the second branch extends from 7 (e'_) to 

7 ( 4 , = t + 1)e . + = 10 s z114 '- 5 


2 rm P c 2 


r tt f 2 F 4 
1 2.3 ( l fc 6,6 


(50) 


(51) 


The largest photon energy measured by LAT, £ ~ 100 GeV, 
corresponds to a pair Lorentz factor 7 ~ 1 0 OO.ZT -11 /F' 2 . 3 , 
which is between j(e'_) and 7 (e(_) (because £_ < 100 GeV < 
£_)_). In further calculations, we make the assumption that the 
LAT power-law spectrum extends well above 100 GeV, imply¬ 
ing that the power-law distribution of pairs extends to 7 > 10 3 , 
as in Ligure 2. A cut-off in the photon spectrum above 100 GeV 
does not affect much the total number of pairs formed but re¬ 
duces the number of higher energy pairs and, consequently, the 
optical flux from pairs formed in-between shocks. 


2.7. Evolution of pair distribution 


The evolution (with observer time t) of the leptons distribu¬ 
tion with energy, Nf)) = dN/dj, is given by 


dU 

~dt 


■N'jiij 




(52) 




£3 

C/3 

*■3 

'3 

Oh 



shock-frame pair energy - y Y 


Fig. 2.— Left panel: distribution with energy (multiplied by 7 ^, to show where the power-law dNi/d'y oc 7 “^ branch occurs) of the pairs formed in 
a relativistic source with same parameters as for Fig 1, calculated numerically by integrating equation (42). Right panel: comparison between the derivative 
d log (dNi / d'y)/d log 7 of the numerical pair distribution and the expected power-law slopes (equations 47 and 48). For e^ (red line), the front is optically thin 

for any photon and the lowest-energy expected branch is not developed. For e & <C all three branches are seen. 


where 

AT raj(7 ) = (53 , 

is the rate at which leptons are created, calculated numerically 
from equation (44) that gives the distribution of leptons formed 
over a dynamical timescale t, and 


d'y dt' d'y 

dt dt dt' 


2r 4g e / 2 
z + 1 3 m e c 


1 ) u 'b(X + 1 ) 


Figure 3 also shows the lepton distribution resulting from 
pair-formation with the above properties and cooled radiatively 
by a magnetic field 

— = 4 (bn) T 2 m p c 2 (56) 

o7T 

that is a fraction b of the post-shock energy. The cooled lepton 
distribution develops a break at an energy 7 C that decreases with 
time, with the lepton distribution being that injected at 7 < 7 C 
and having a slope larger by 1 than that injected for 7 > 7 c . 


= -M7 2 - 1) , k r = 1 *L.TB 2 (Y + 1) (54) 

dir(z + 1 ) m e c 

is the radiative cooling rate of pairs, t' being the shock-frame 
time, u' B = B 2 /8tt the shock-frame magnetic energy density, 
and 


d f‘°° rr r°° 

Y =31 ^“=3^/ (55) 

is the Compton parameter (the ratio of the inverse-Compton 
to synchrotron losses), r sc = cr e V/( 47 ri? 2 ) being the optical 
thickness to photon scattering by N leptons in a source of ra¬ 
dius R. 

The lepton distribution Mj ;n y depends on the properties of 
the high-energy emission (<f>, £>,, a, 3) and of the photon source 
(T, R). Figure 3 shows the instantaneous and integrated in¬ 
jected lepton distributions Af. in j for a source with constant 
a, (3, a source deceleration corresponding to a blast-wave in¬ 
teracting with a wind-like medium - F = r o (t/i 0 ) -1 / 4 , and 
an evolution of the high-energy spectrum break as expected for 
the forward-shock emission: £5 = £b(t 0 )(t/t 0 )~ 3 ^ 2 . Because 
£b < (equation 19), the photon front is optically-thick to 
pair-formation above £_ (equation 21 ), thus the pairs have a 
Afi n j oc 7 -fi distribution above 7 (e(_) (equation 48), which in¬ 
creases with time: 7 (e , _) oc F 5 f 2 oc f 3 / 4 . 


That feature for the radiative cooling of a power-law distri¬ 
bution of particles can be derived from the kinetic equation for 
particle cooling (equation 52), rewritten as 


dt 

using equation (54). 
a (t) 7 -p , leads to 


= fci7 + k, 


9(fA0 

d'y 


(57) 


Trying a power-law solution Af = 


0 = P + ( p _ 2 )^a7 1 v (58) 

dt 


If the first term on the right-hand side is dominant (i.e. 
da/dt krd'y and the radiative cooling term is negligible), 
then p = (3 and da/dt = ki, thus a = J kidt and Af = 
(f k i dt) r y ~^, hence the effective distribution is the integrated 
injected distribution, at energies 7 that satisfy 7 <C 7 cL where 


7 


(i) = 


da/dt 

k r a 


h 

kr f kidt 


(59) 


If the second, radiative cooling term on the rhs is dominant, 
then p = ft + 1 and a ~ ki/k r , hence Af = (ki/k r ) 7 ~(/ 3 + 1 ). 
This solution exists for \da/dt\ k r a 7 , which is equivalent to 

7 7 c, where 




\da/dt\ 

k r a 


1 d ki 


ki 


dt k r 


( 60 ) 


























pair LF - y 
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Fig. 3.— Solid lines show the distribution of leptons formed at three epochs (10 s, 100 s, 1 ks - colored lines) and total distribution at 2 ks (black), integrated over 
the indicated source evolution and without radiative cooling (b = 0). The quasi-instantaneous distributions are normalized to that integrated up to 2 ks. The source 
GeV fluence is constant, its deceleration T(t) and break energy e^ evolutions correspond to a blast-wave interacting with a wind-like external medium. Black lines 
show time-integrated (up to 2 ks) lepton distribution with synchrotron and inverse-Compton cooling, for different fractions b of the post-shock energy that is in the 
magnetic field. The cooled lepton distributions display a break at the energy for which leptons lose half of their energy over a dynamical timescale. Above the break, 
the lepton distribution slope increases by unity (right panel). The higher the magnetic field parameter b, the lower the lepton cooling-break energy. 


The above two cooling energies are comparable 

7c 1 ) _ ht ! din \kj/k r \ 

-y( 2 ) J kjdt / dint 

if ki does not vary too fast and if ki/k r is moderately evolving 
in time. Thus, the effective lepton distribution is 


•A/" ( 7 ) = 


(/ hdt)"f P 7 < 7 ^ X) 

|i 7 -(/3+i) 7 

li/'p 


(62) 


For a constant injection rate ki, the 7 c becomes simpler 


1 


7c 


1 


k r t TB 2 {Y + l)t 


oc 


(63) 


and the effective distribution can be approximated as 


with 7 the energy of a lepton that had initially an energy j 0 , 
after a timestep St. Pair-energy tracking means accounting for 
that 

(1) a fraction min[l, 5t/t c (^)\ of the SN = Af^S'y leptons 
existing in a cell (7,7 + £7) exit that cell due to their cooling 
during St, 

(2) a fraction vnox[\,t c { r )) / St] of the leptons SNi n j = 
Ninj{'y)S'ySt injected in a cell remain in that cell after cool¬ 
ing for St, 

(3) leptons existing in all cells between energies 70 ( 7 ) and 
7 o(7 + £7) cool during St to cell (7,7 + Sj), 

(4) afraction [t c ( 7 0 —> 7 ) —f c (7 o —> 'y + Sy)]/St ofthe leptons 
injected at energy 7 0 cool to cell ( 7,7 + S'y) during St, with 
f c ( 7 o —i► 7 ) = ( 7 -1 — 7 tT 1 )/^c. the cooling time from energy 
7 o to 7 . 


3. RADIATION EMISSION 


•A/"(7) - A/"(7c) 



7 < 7c 
7c < 7 


N (7c) = y 1 


-03+1) 


(64) 


We note that the cooled lepton distribution is not calculated 
from equation (52) because it is unstable and ’’suffers” from 
a Courant-like condition, with the timestep St being upper- 
limited by the cooling time t c ( 7 ) = (k r r ))~ 1 of the highest 
energy leptons 7 max in the calculation, which can be t c <C t 
for bn-ymax > 10 3 r^ 3 (y + 1) _1 . Instead, the lepton distribu¬ 
tion is calculated numerically by tracking the flow of particles 
on a 1 -dimensional energy grid, using the cooling law of equa¬ 
tion (54) 


7 + 1 
7-1 


lSL±l e ^rSt > I 7^1 1 + krSt 

7o - 1 7 7o 


(65) 


The calculation of the synchrotron self-Compton from 
pairs three components: synchrotron emission, synchrotron 
self-absorption, and inverse-Compton scattering of the self- 
absorbed synchrotron spectrum. We consider only the first 
inverse-Compton scattering, which is appropriate approxima¬ 
tion when the Compton parameter Y is sub-unity and a nec¬ 
essary approximation when leptons radiating at the observing 
frequency v cool faster than they are created. The former case 
requires a magnetic field in the shock fluid that is not much be¬ 
low equipartition (b > 1 (K :i ); for the latter case, the argument 
is that, in an source that is optically-thin to electron scattering, 
the lepton distribution can change substantially during the time 
it takes a photon to cross the source and be scattered, hence 
the lepton distribution that produced the seed photon is not the 
same as the lepton distribution that upscatters it, an a time- 
dependent treatment of upscatterings is needed. For Y > 1, 
ignoring higher-order scatterings leads to an underestimation 
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of the flux at higher energies (above X-rays) and an overesti¬ 
mation of the synchrotron and first inverse-Compton flux above 
the cooling frequency (which would be lower if higher-order 
scatterings were accounted). 


3.1. Synchrotron Emission 


For a relativistic source moving at Lorentz factor F, the re¬ 
ceived synchrotron flux from a distribution of leptons Af( 7) at 
observer frequency v is 


F sy(y) = J d l^{l) P 'sy ( 66 ) 


where the relativistic boost of the comoving-frame emission at 
frequency v' = (z + l)jz/r gets only one power of T from 
the contraction of photon arrival time dt = dti a b/T 2 = dt'/T 
relative to the comoving-frame emission time dt' (the boost F 
in photon energy is ’’lost” to that a comoving energy range dv' 
is stretched into F dv' for the observer, and the angular beam¬ 
ing boost r 2 is ’’lost” because, for the observer, that beaming 
reduces the solid angle of a spherical source by a factor I’ 2 ), 
and 


KyV, 7) 


e 3 B 

m e c 2 



(67) 


is the specific synchrotron power for a lepton, with f sy the ’’syn¬ 
chrotron function” and 


u 'sy(n) = 7R — B 7 2 = 3-3 x lo6 £?o7 2 Hz ( 68 ) 
lb m e c 


satisfies v sy ( 7 ^) = v. Owing to the exponential cut-off of the 
synchrotron function at x > 1, only leptons with energy above 
7 „ produce the synchrotron emission at frequency v. Then, 
approximating the synchrotron function by only its 2 : < 1/2 
branch, one obtains 


F S y{l') 


0.92 x IO - 57 

(Z + l) 8 / 3 



dj 


TV/t) 

y 2/3 


Jy 

(73) 


For a power-law distribution of particles, /V(7) oc 7 p , this 
leads to the well-known spectrum F sy (y) oc v~( p ~ lp2 . 


3.2. Synchrotron Self-Absorption 


Starting from equation (6.49) of Rybicki & Lightman (1949), 
taking into account that the lepton column-density is N/4ttR 2 , 
and using the synchrotron emissivity per lepton given in equa¬ 
tion (67), the synchrotron self-absorption optical thickness at 
observer frequency v = To’ /(z + 1) is 

Ta ^ = ~32 t t 2 R 2 


e 6 B 


(Z + l)v\ ilTl e c) 


(74) 


X 



d m 7 )\ 

dl \ 7 2 ) 


Quick progress toward a simpler form can be made if we retain 
only the x < 1/2 branch of the synchrotron function (equation 
70) and set f sy (x > 1/2) = 0 (motivated by the exponential 
cut-off): 


is the synchrotron characteristic frequency at which a lepton of 
energy 7 radiates most of its emission. The synchrotron func¬ 
tion is an integral over the modified Bessel function of 5/3 order 
and has the following asymptotic behavior 


r 1.71 x i > 3 £<1 

\ \.2'Sx x ! 2 e~ x x 1 


(69) 


This approximation would be useful if the synchrotron emis¬ 
sion at frequency v’ were produced by leptons whose charac¬ 
teristic synchrotron frequency 1 / is far from however, the 
opposite is true. We approximate the synchrotron function with 
the asymptotic behaviors given above but with coefficients such 
that i) f sy is continuous at x = 1/2 (where x 1 ^ 2 e~ x has a 
maximum) and ii) its integral is equal to that of the exact syn¬ 
chrotron function, J f sy {x)dx = (4/3) 3 , yielding a power-per- 
lepton P'(j) = J P' sy (y', f)dv' = (4/3)a e c^ 2 (B 2 /8tt). The 
following approximation 


f 1.50 x 1 / 3 x < 1/2 

\ 2.77x 1 ^e~ x x > 1/2 


(70) 


satisfies the above constraints and has a maximum f sy ( 0.5) = 
1.2 (that of the exact synchrotron function is f sy { 0.3) = 0.92). 

Substituting equations (67) and (68) in (66), we get 


4.3 xlO- 56 
SS() " (3 + 1) 3 

in cgs units, where 


dlAf{l)fsy (|) Jy (71) 


y v = 5.5 x 10 


,-4 


(z + l)l 

TB 


11/2 


(72) 


Ta{v) 


> 


0.0083 

R 2 


-,5/3 


(z + l)^m e c 


(e 4 B) 


2/3 



The integral above can be re-written as 


4 7 


(z + l)v. 


5/3 


B 2/3 


(75) 


(76) 



in cgs units. For a power-law distribution of leptons with 
energy, Af( 7 ) oc y~ v , the ratio of the two terms above is 
(3 p + 2)/4, where p = /?,/?+ 1 (equation 62 and Figure 3). 
For a typical spectrum of the high-energy photons, /3 ~ 2, this 
ratio is 2 or 11/4, thus the first term is dominant. Either term is 
proportional to 7,7^ p+2 ^ ) and yields r a (y) oc ^ - (p+ 4 )/ 2 . 

The emerging synchrotron flux is the intrinsic flux (given in 
equation 71) reduced by the escaping fraction (equation 24): 


F^ bs \u) 


1 - e ~ Ta ^ 
4«W 


Fsyiy) 


(77) 


Below the self-absorption frequency v a , where T a {y a ) = 1, the 
received synchrotron flux satisfies Fs y bs \v) = F sy (v)/T a (v) 
which, for a power-law distribution of particles, is Fg y bs ^ {v < 
v a ) oc i/-(p- 1 )/ 2 / J y-(p+ 4 )/ 2 = 1 / 5 / 2 5 another well-known re¬ 
sult. 
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3.3. Inverse-Compton Emission 


For a lepton of energy 7 scattering a photon of frequency 
v' Q in the Thomson regime, the inverse-Compton emissivity at 
photon frequency v' is (equation 7.26a in Rybicki & Lightman 
1979) 


Uv') 


3 CJeF' 0 ( v' 

4 (7^r c V 4 tM 


(78) 


where f )' is the comoving-frame energy flux of;/ photons and 


, , x f 2a; In x — 2a; 2 + a: + 1 (I67 4 ) 1 < x < 1 

= { 0 otherwise 

(79) 

By integrating ■N’{y){i , ' 0 )j'. ic (i'') over the lepton distribution and 
over the incident synchrotron photon spectrum, one obtains 
the comoving-frame inverse-Compton luminosity, which is en¬ 
hanced by a factor T to yield the lab-frame inverse-Compton 
luminosity. Then, the received inverse-Compton emission from 
pairs at frequency v = Tv '/ (z + 1) is 


F ic {v) 


3(2: + 1) 

8ndf 


oylV 


d'yNipi) 





The comoving-frame flux of incident photons F' 0 {y' 0 ) in a 
source of radius R is related to the received synchrotron flux 
at photon frequency v 0 = Tv' 0 /(z + 1 ) through 


Fsy(vo) = Z -^T{A-kR 2 )F' 0 {v' 0 ) (81) 


Substituting in equation (80), and changing from comoving- 
frame to observer-frame photon frequencies, it follows that the 
received inverse-Compton flux can be calculated from the re¬ 
ceived synchrotron flux (equation 71): 


where r e = (o e Ni)/(4nR 2 ) is the pairs’ optical-thickness 
to photon scattering. Equation (84) means that the inverse- 
Compton specific flux is the integral over the lepton distribution 
of the scattered synchrotron flux. 


4. EMISSION FROM PAIRS FORMED IN THE SHOCKED FLUID 


4.1. Approximate dependences 


The most important parameters that determine the pair emis¬ 
sion are those that set the number of pairs - the blast-wave ini¬ 
tial Lorentz factor r o and the afterglow high-energy fluence $ 
- and the magnetic field - the nb product (equation 56 ). Less 
effective, but still relevant, are three other parameters that deter¬ 
mine the number of pairs: the slopes a and 7 of the high-energy 
spectrum, and its break energy Si,. 

Equations (71) and (84) suggest the following dependences 
for the synchrotron and inverse-Compton flux from pairs: 

N+ Nl 

F sy (v) oc BTN± , F ic {v) oc —F sy oc -^ST ( 85 ) 


where R oc T 2 f (equation 6 ), B oc VnbY (equation 56), and 
N± oc $ 2 /(r 6 f 2 ) (equations 30 and 34). Equation (71) ac¬ 
tually means that F sy (y) oc N(> 7^), hence the use of N± 
here is accurate only when the number of pairs above 7 ,, is a 
fixed fraction of the total number of pairs. That is satisfied only 
above the pair cooling-energy, where pairs produced during one 
cooling timescale reside. With the above substitutions, we find 


that 


F sy (v ) oc 


^(nh) 1 / 2 

TH 2 


F ic {y) oc 


^(nb) 1 ^ 2 

r 14 t 6 


( 86 ) 


For a wind-like medium, where n oc R 2 and F oc t 7 4 , we 
arrive at 


Fic(v) 


3o e 

167tR 2 



where Fj y bs •* is the self-absorbed synchrotron flux (equation 
77), because synchrotron self-absorption reduces the flux of 
photons incident on a scattering lepton in the same way that 
it affects the received synchrotron flux. 


An approximation that changes very little the spectrum of 
the inverse-Compton pair emission is that where photons of fre¬ 
quency v result from the upscattering by a lepton of energy 7 of 
only synchrotron photons of frequency v a = 3 h/( 47 2 ), which 
is motivated by that the average energy of an upscattered photon 
is ( 4 / 3 ) 7 2 u 0 . This is equivalent to approximating the inverse- 
Compton function with a (5-function 
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Then, equation (82) becomes 
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This shows that the inverse-Compton flux has a very strong 
dependence on the high-energy afterglow fluence and a super¬ 
strong dependence on the Lorentz factor of the GeV afterglow, 
which suggest that inverse-Compton emission from pairs could 
be relevant (i.e could overshine the synchrotron flux) only for 
the brightest LAT afterglows (<!> ~ 10 - 5 ergcm~ 2 ) and for 
the slowest GeV sources in which pairs are still optically-thin 
(F x) 130 - equation 32). 


4.2. Light-curves and spectra 

The pair light-curves shown in Figure 4 illustrate the corre¬ 
lation of the pair flux with the observable LAT fluence $ and 
the unknown source Lorentz factor F 0 at the peak epoch t„ of 
the LAT light-curve. Those light-curves were obtained by in¬ 
tegrating the synchrotron and inverse-Compton fluxes given in 
equations (71), (76), (77), and (82) (or 84), over the decelera¬ 
tion of a blast-wave interacting with a massive-star wind. Al¬ 
though N± oc 4> 2 /{T 6 0 t 3 0 ) is satisfied by the numerical calcu¬ 
lation of the pairs formed, the numerical pair fluxes display a 
weaker correlation with $ and F„ (and also with t Q ) than given 
in equation (87), which is due to the use of N± in the derivation 
of that equation. 
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Fig. 4.— Dependence of pair emission on the fluence of the high-energy LAT afterglow fluence <f> (left panel) and on the GeV initial Lorentz factor r o (right 
panel), with other parameters as indicated. Red lines are for synchrotron and blue lines for inverse-Compton. 


Equations (30) and (34) show that the number of pairs is 
weakly dependent on the unknown break energy £& of the LAT 
spectrum (another model parameter ), with more pairs being 
formed for a lower £&, because that increases the optical thick¬ 
ness to pair formation (equations 15 and 21). The brightness of 
the LAT high-energy spectral component at sub-MeV photon 
energies relative to that of the burst is the criterion for choosing 
the two prescriptions given in Figure 5 for the unknown break- 
energy: GRBs with a fast-decaying tail require that the £& of a 
bright LAT component remains above MeV for the duration of 
the tail, while bursts with a slowly-decaying tail allow lower £& 
(decreasing or not). As expected, a lower £/, yields a brighter 
pair emission, and a decreasing £5 leads to a slower decay of 
the pair light-curve. The latter behavior provides a criterion for 
identifying early optical afterglows produced by pairs: slowly- 
dimming pair afterglows (due to a decreasing £ 5 ) cannot follow 
fast-falling bursts (which are incompatible with a decreasing 
£ 5 ). However, fast-falling pair afterglows can follow either type 
of burst tail (fast or slowly decreasing). 

4.3. Application to GRB 130427A 

Figure 6 shows a fit to the super-bright optical flash of GRB 
130427A (RQD2/Raptor - Vestrand et al 2014) with the syn¬ 
chrotron emission from internal pairs formed from the high- 
energy emission monitored by LAT (Fan et al 2013, Tam et al 
2013, Ackermann et al 2014). Observations set the high-energy 
fluence $ and the spectral slope /3 above the unknown break 
energy £&, which is a free model parameter. The initial Lorentz 
factor T 0 of the high-energy photons source and the magnetic 
field parameter b in the shocked fluid are two other model pa¬ 
rameters. 

As indicated by equation (87), with the high-energy flu¬ 
ence 4) set by observations, the brightness of the optical flash 
of 130427A constrains the combination b/Tl 2 . We find that 
T 0 < 300 is required to match the brightness of the 130427A 
optical flash because, for higher Lorentz factors, the number of 


pairs formed is too small, and the maximal optical flux from 
pairs, obtained for a magnetic field that brings the peak of the 
self-absorbed synchrotron spectrum in the optical, falls short of 
the peak brightness of GRB 130427A’s optical counterpart. 

The decay of the pair synchrotron optical light-curve depends 
primarily on the LAT light-curve 4>(t) and on the blast-wave de¬ 
celeration T(f). As those quantities are already set, the optical 
flash decay (F 0 oc t~ 2 ) constrains the slope of the LAT spectral 
component below its £& break and the evolution of £ 5 . 

For a = 2/3 (i.e. £t is the peak energy of a synchrotron 
or inverse-Compton spectrum without cooling), we find that 
£b oc t [/ “ 2 is required to match the optical light-curve decay 
at 10-100 s. This time-dependence is consistent with the be¬ 
havior of the cooling-break of the synchrotron spectrum from 
the forward-shock (and for a wind-like medium), but a = 2/3 
is inconsistent with the expected value a = /? — 1/2 in that 
case. 

For a = 3/2 (when £5 would be the injection peak of the 
sy/ic spectrum with electron cooling) or for a = j3— 1/2 = 1.7 
(when £5 would be the cooling-break of a sy/ic spectrum), we 
find that the optical flash decay requires £& oc f 2 , which is con¬ 
sistent with the evolution of the cooling-break of the inverse- 
Compton spectrum from the forward-shock (and for a wind-like 
medium). 

Therefore, fits to the decay of the prompt optical emission 
of GRB 130427A with emission from internal-pairs sets con¬ 
straints on the unmeasured peak energy £b of the LAT spectral 
component that do not elucidate its shock origin. Furthermore, 
numerical fits to the multiwavelength emission of this afterglow 
show comparable contributions to the LAT emission arising 
from both synchrotron reverse and forward shocks (Panaitescu 
et al 2013). 

Figure 6 also shows that the X-ray emission from pairs and 
that LAT component contribution to the X-ray are below the 
fluxes measured by Swift, and that the formation of enough 









iMill iiiiiiij iiniiij iiniiij iiiiiiij iiiiiiij niiiiij iiiiiiij ma 


e b = 10 5 (t/10)" 3/2 eV 


1 10 


5 10 


l 10" 


o 

o 



time (s) 


10 10 10 10 10 
frequency V (Hz) 


Fig. 5.— Left panel: optical light-curves from pairs formed in the shocked fluid, for the indicated parameters, and for the wind-like external medium given in 
equation (35). Solid lines are for a constant break-energy e & of the high-energy emission that forms pairs. The synchrotron flux (red line) exhibits a decay that is 
slightly steeper than estimated from equation (87), while the inverse-Compton flux (blue line) has a decay that is significantly slower. Dashed lines are for a decreasing 
£&, which yields a slower flux decay. Right panel: synchrotron and inverse-Compton spectra at t = 100 s, for a decreasing Dotted lines show spectra without 
accounting for synchrotron self-absorption (SSA), solid lines are for absorbed spectra. Owing to the large magnetic field parameter b, the Compton parameter is below 
unity and the inverse-Compton emission is dimmer than synchrotron at all frequencies of interest. The synchrotron cooling break is slightly below the optical. 


pairs to produce a bright optical flash does not entail a high 
attenuation of the LAT spectrum above 10 GeV. For the high¬ 
est Lorentz factor r o that allows a good fit to the optical flash, 
the intrinsic power-law spectrum above 10 GeV is attenuated 
by up to 50 percent at t = 10 s (when attenuation is maximal), 
which is not inconsistent with the detection by LAT of a 70 
GeV photon at 18 s. However, Lorentz factors F 0 < 200 are 
incompatible with that detection. 

Although a good fit with the internal-pair emission for the 
optical flash of GRB 130427A is obtained, we do not propose 
this origin for the optical counterpart of GRB 130427A because 
modeling of the broadband (radio, optical. X-ray, and GeV) 
emission of this afterglow (Panaitescu, Vestrand & Wozniak 
2013), from 100 s to tens of days, has shown that its wind-like 
ambient medium must be very tenuous, which leads to an initial 
Lorentz factor L„ ~ 750 that is much higher than allowed by 
fitting the optical flash with internal-pairs emission (r o < 300). 

Synchrotron emission from external pairs formed ahead of 
the blast-wave from burst MeV photons scattered by the am¬ 
bient medium and, then, accelerated by the forward-shock can 
also produce a bright optical flash (R < 10) lasting for 100 s, 
provided that the initial source Lorentz factor is l’„ ~ 200 (fig¬ 
ures 4 and 7 of Kumar & Panaitescu 2004). Vurm, Hascoet & 
Beloborodov (2014) have found that the optical flash of GRB 
130427A can be explained with synchrotron emission from ex¬ 
ternal pairs accelerated by the forward-shock if that shock’s 
Lorentz factor is a low L = 200. 

A similar model, but not investigated here, is the emission 
from the shock-accelerated external pairs formed from after¬ 
glow MeV-TeV photons ahead of the forward-shock. In one 
variant of that model - pair-formation from unscattered GeV 
photons - the number of pairs is strongly decreasing with the 
source Lorentz factor, therefore it requires a low r o < 300 to 


account for the optical flash of GRB 130427A. In the other vari¬ 
ant - pair-formation from GeV photons scattered by the ambi¬ 
ent medium (which decollimates photons sufficiently to lower 
significantly the pair-formation threshold-energy and and en¬ 
riches with pairs the medium ahead of the blast-wave) - the 
number of pairs should be less dependent on r o . Owing to its 
similarity to the pair-wind formed from scattered burst MeV 
photons, this model may also require a low r o to account for 
the optical flash of GRB 130427A. 

If all pair-based models for this flash require low Lorentz 
factors (for the seed-photon source) that are incompatible with 
the afterglow T inferred from multiwavelength data modeling, 
the bright optical flash of GRB 130427A should be attributed 
to the reverse-shock (Meszaros & Rees 1997) that energizes 
some incoming ejecta in an initial injection episode, followed 
by a quiet period when the ejecta electrons cool radiatively and 
yield a fast-decaying flux, followed by a second, longer-lived 
injection episode, during which the reverse-shock produces the 
optical emission measured for the early (up to few ks) afterglow 
of GRB 130427A (as proposed by Vestrand et al 2014). 

5. CONCLUSIONS 

In GRB afterglows, test photons of lab-frame energy above 
~ 10 MeV form pairs in interactions with target photons that 
are above the threshold for pair-formation. The number of 
pairs depends moderately on the unknown break-energy of 
the high-energy component (LAT measures only photons above 
£b), strongly on the afterglow GeV output (which is the observ¬ 
able LAT fluence $), and very strongly on the Lorentz factor T 
of the GeV source. 

Below the radiative cooling break, the brightness of the syn¬ 
chrotron emission from (internal) pairs, formed in the shocked 
fluid (between the reverse and forward shocks), depends on 
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Fig. 6. — Left panel: early optical and GeV emission for GRB 130427A and a fit to the optical flash (up to 100 s) with synchrotron emission from internal pairs 
formed in the shocked fluid. LAT fiuence (model input) after the t 0 = 10 s peak is indicated; LAT spectral slope is /3 = 2.2 (Tam et al 2013). Model parameters 
are: source Lorentz factor before LAT peak r o = 300, break-energy of the LAT spectrum = 500 (t/t 0 ) 1 /2 keV , magnetic field parameter b = 10“ 3 . The burst 
ambient medium is that of equation (35), which sets the source dynamics: T(f) = r o (t/t 0 ) — 1 ' 4 . After 100 s, the optical emission is well-fit by the reverse-shock 
emission (Panaitescu et al 2013). Mid panel: for the parameters of the optical flash fit, the pair emission is dimmer than the early X-ray emission of this burst, 
monitored by Swift’s BAT and XRT. The 10 keV emission from the LAT component with an assumed low-energy slope a = 2/3 below St, is also shown. Part of the 
burst tail and early X-ray afterglow can be explained with the forward-shock emission. Right panel: photon-photon attenuated spectra (dashed lines) at the LAT peak 
light-curve epoch, when attenuation is maximal, for the measured LAT spectrum (solid line) and various initial F 0 . For r o = 300, a moderate absorption occurs 
above 10 GeV, corresponding to a flux reduction of at most 50 percent. For smaller r o , attenuation is stronger. 


their number (set by one observable - $ - and two model pa¬ 
rameters - r and £b) and on the strength of the magnetic field 
between shocks (a third model parameter). For an intermedi¬ 
ate/low r, pairs produce bright optical early afterglows even 
for a magnetic field that is several orders of magnitude below 
equipartition. In fact, strong magnetic fields do not warrant a 
much brighter optical emission because an enhanced radiative 
cooling reduces the number of pairs of sufficiently high energy 
to radiate synchrotron emission in the optical. 

The correlation between the number of pairs and the attenu¬ 
ation of the LAT spectrum, induced by the dependence of these 
two features on T, provides a way to identify optical coun¬ 
terparts that originate from internal pairs (formed in the GeV 
source). For the most relativistic afterglows (F > 500), the 
internal-pairs emission should be dim and the LAT spectrum 
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should be an unattenuated power-law, both because few pairs 
are formed. 300 < T < 500 yields a moderately bright op¬ 
tical flash and no detectable attenuation of the LAT spectrum. 
For the less relativistic afterglows (100 < T < 200), when 
many pairs are formed, there should be a bright optical emis¬ 
sion from pairs, accompanied by a significant attenuation of the 
LAT spectrum above 1 GeV. 

An additional criterion for identifying optical counterparts 
from internal-pairs emission arises from that GRBs with fast- 
decaying tails require a peak energy £ b > 10 MeV of the 
LAT spectrum, which yields dimmer and faster-decaying op¬ 
tical emission from pairs. Slowly-decaying GRB tails do not 
exclude bright optical flash from pairs, hence there should be 
some correlation between the speed of the GRB tail decay and 
the brightness of the pair optical flash. 
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